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Abstract 

Besides being an ancient pigment, BaCuSi206 is a quasi-2D magnetic insulator with a gapped spin 
dimer ground state. The application of strong magnetic fields closes this gap creating a gas of bosonic spin 
triplet excitations called triplons. The topology of the spin lattice makes BaCuSi206 an ideal candidate 
for studying the Bose-Einstein condensation of triplons as a function of the external magnetic field, which 
acts as a chemical potential. In agreement with quantum Monte Carlo numerical simulation, we observe a 
distinct lambda-anomaly in the specific heat together with a maximum in the magnetic susceptibility upon 
cooling down to liquid Helium temperatures. 

PACS numbers: 75.45.+j, 75.40.Cx, 05.30.Jp, 67.40.Db, 75.10.Jm 
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About 2000 years ago, early Chinese chemists synthesized barium copper silicates for the first 
time, and used them as pigments forpottery and trading as well as for large Empire projects 
such as the Terracotta Warriors |lll|2l|3|], preceding even the invention of paper and the compass. 
BaCuSi206, also known as Han Purple |ll|], is then possibly the first man-made compound con- 
taining a metallic bond. In its layered crystallographic structure, pairs of Cu^+ ions form dimers 
arranged in a square lattice, with Cu-Cu bonds projected normal to the planes Neighboring 
Cu bilayers are weakly coupled, i.e. the magnetic system is quasi-two dimensional Applied 
magnetic fields in excess of 23. 5T suppress the spin singlet ground state in BaCuSi206, giving 
rise to a gas of quantum spin-triplet excitations. 

At low temperatures, an ideal gas of bosons undergoes a phase transition into a condensate 
with macroscopic occupation of the single-particle ground state. A great deal of interest in this 
phenomenon was triggered by the discovery of the anomalous behavior of liquid helium. Upon 
cooling, liquid ^He exhibits a A-transition in the specific heat at 2.17 K that signals the onset of a 
zero-viscosity superfluid state 1^. Superfluidity arises from the macroscopic fraction of He atoms 
that occupy the single particle ground state. Within the last decade, Bose-Einstein condensation 
(BEC) was also realized for dilute clouds of atoms at temperatures lower than ten millionths of a 
degree Kelvin The possibility of producing BEC with quantum spin magnets [8] has stimulated 
considerable experimental effort to find the candidate materials ||9[l. 

Spin-dimer systems such as SrCu2(B03)2 and TlCuCla 111 111 have recently become of in- 
terest owing to the ability of strong magnetic fields to generate a gas of = 1 spin-triplet states 
moving in a non-magnetic background. These triplons can be regarded as bosonic particles 
with hard core repulsions that carry a magnetic moment, but no mass or charge. The external field 
plays the role of a chemical potential in controlling the number of particles. The Ising compo- 
nent of the inter-dimer exchange interaction generates an effective repulsion between neighboring 
bosons. 

Upon cooling, the gas of triplons can either crystallize or condense in a liquid state, depending 
on the balance between the kinetic energy and the repulsive interactions [9]. In particular, if 
the kinetic energy dominates and the number of triplons (or total magnetization) is a conserved 
quantity, the system undergoes BEC corresponding to the coherent superposition of S'^ = 1 
spin-triplet and S = sin glet states on each and every dimer. The BEC of triplons was proposed 
to occur in TlCuCk Hill Q and KCuCls fly. However, recent data suggest that the 
anisotropic crystal structure of these spin-ladder systems breaks the U(l) rotational invariance that 
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is required to have gapless Goldstone modes Hence, the Goldstone modes acquire a finite gap 
or mass. In this report we introduce a new system whose bi-layer structure provides a reahzation 
of a quasi-two dimensional gas of triplons. 

The range of fields and temperatures required to obtain a BEC of triplets in BaCuSi206 creates 
a unique situation from the experimental standpoint: the triplons dominate the thermodynamics 
in a very simple and predictable fashion. The simplicity of BaCuSi206 is best described by a 
spin Hamiltonian that includes only the nearest-neighbor Heisenberg antiferromagnetic exchange 
couplings. The Cu^+ ions and the Si04 tetrahedra in BaCuSi206 are arranged in layers parallel to 
the (001) crystallographic plane, as shown in Fig. 1(a). Within each Si-O-Cu layer, the Cu^+ ions 
form a square lattice bi-layer of = 1/2 spins The relevant Heisenberg exchange couplings 
within the bi-layers are the intra-dimer interaction, J, and an inter-dimer nearest-neighbor inter- 
action J' (Fig. 1(a)). The different bi-layers are antiferromagnetically coupled via the effective 
exchange constant The resulting spin Hamiltonian is: 



Hg — J S[ i • Si^2 + Si^2 • Si+ea,2 



i,OL,f3 



+ J" ^ Si^i • Si+^,2 - Qw^^bH ^ Slf^ 

i i,/3 



(1) 



where a = {1, 2} is the direction index, [3 = {1, 2} is the layer index, ei = x, 62 = y and z 
are unit vectors along the crystallographic axes, and g\\ = 2.306 ± 0.008 is the component of the 
gyromagnetic tensor along the c-axis, which is the direction of the applied magnetic field H. The 



value of J = 4.5meV is known ff 
J'' < J' <C J. Precise values of 




studies. [5], which also show that 
he experiments described below. 



FIG. 1: (color online) a) Cu^+-dimer plane in BaCuSi206 in a square lattice arrangement. J and J' are the 
intra-dimer and inter-dimer exchange constants, b) Evolution of the isolated dimer spin singlet and triplet 
with an applied magnetic field H. c) Partially occupied hard core boson lattice. Occupation (chemical 
potential) is given by the external magnetic field, d) Transmission-light picture of a BaCuSi206 single 
crystal. 
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In the limit J' = J" = 0, the ground state of is a product of local singlet-dimer states: 
|*o) = where 10?) is the singlet state between two spins on a dimer i. J = 4.5 meV 

is the energy gap for a spin triplet excitation in a single isolated dimer This gap decreases 
linearly with an applied magnetic field as shown in Figure 1(b). Since J' <C J, the spectrum 
is dominated by the two low energy states of the isolated dimer (the = singlet and the 

= 1 triplet) once gufisH becomes of the order of J, By projecting H onto this low energy 
subspace, we derive an effective Hamiltonian, Hef, in which each dimer is represented by one 
effective site with two possible states. In this way, we are neglecting the effect of the high energy 
triplet states and, therefore, the highest order terms of Heff are linear in f and = (first order 
perturbation theory). A similar procedure was used before by Mila [18] to describe spin ladders 
in a magnetic field. We associate the two low energy states with the two possible states of a hard 
core boson on a lattice ||5|] (see Fig. 1(c)). The empty site corresponds to the singlet state \(pl), 
while the site occupied by the hard core boson represents the = 1 triplet state: 10?) |0)i, 
|(/)|, = 1) ^ 6j |0)i,where |0)i represents the empty state at the site i. In this language, the 
effective low energy Hamiltonian is: 



where the chemical potential is /x = J — gwii^H, t = V = J' /2, and t' = V = J" j^, Eq.|2| 
describes a gas of hard core bosons with nearest-neighbor hopping and repulsive interactions. For 
a less than a critical field i/ci {g\\^B^c\ = J — '^J^ — J^' /2), the ground state of Hejf is completely 
empty, i.e. all dimers are spin singlets, because the chemical potential, /x, is large and positive. A 
finite concentration of particles p = m (m is the magnetization per site) emerges in the ground 
state only when H exceeds Hd, Because the condition V < 2tis not fulfilled, the triplets or hard 
core bosons never crystallize for any concentration p < 1, 

The single crystal samples used in this study were prepared using the floating-zone technique. 
BaCOa, CuO, and Si02 were used as starting materials and the powders of stoichiometric compo- 
sitions were mixed and sintered at 700° C for 15h. The sintered powder was reground and molded 
into a rod under a pressure of 400kg/cm^ for about lOmin. The rod is sintered between 1000 and 

n 

1100°C for about 100 h. The sintered rod is then used for floating zone method II22I1 . The floating 
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FIG. 2: (color online) a) Magnetization vs field at different temperatures from 0.4K to 40 K, measured in 
both a 65 T short-pulse, and a 45 T mid-pulse, capacitor driven pulsed magnets, b): Magnetic susceptibility 
vs temperature from measurements in A, and from our model using a quantum Monte Carlo algorithm. 

zone grov^th was done tv^ice. In the first run the speed of the grovv^th was set to 50mm/h, and in the 
second run to ~ 0.5 mm/h, both in O2 flow of ^ 200 cc/min. The purple color and texture of our 
single crystal sample is appreciated in a picture taken in transmission (backlight) mode (Fig. 1(d)). 
Figure 2a show^s the magnetization versus field data obtained in capacitor-driven pulsed magnetic 
fields The magnetization increases continuously betv^een the critical field Hd = 23.5 T and 
the saturation field Hc2 = 49 T at T = 0.4 K. The nearly linear slope of the magnetization contrasts 
sharply the staircase of plateaux observed in SrCu2(B03)2, where the gas of triplets is known to 
crystallize in a structure that is commensurate with the underlying crystallographic lattice 
Pulsed magnetic field experiments on BaCuSi206 were repeated at different temperatures up to 
40 K and cross-calibrated against low-field SQUID magnetometer data The structure ob- 
served around 48 T at low temperatures, which we could not reproduce in later experiments, is 
likely an artifact caused by mechanical vibration of the probe. A plot of the slope x = dM/dH 
of M versus H (magnetic susceptibility) at i7 = 37 T is shown in figure 2(b) together with the 
results of a Monte Carlo simulation of JLff in a finite lattice of size L^, with L = 12, performed 

nrt 

using the directed-loop algorithm ||20l|21I] with parameters: J = 4.45 meV, J' = 0.58 meV and 
J'^ = 0.2 = 0.116 meV. The agreement between experimental data and model calculations is 
remarkable. Notably, the rapid drop of the magnetic susceptibility below 4 K indicates the onset 
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of an ordered low temperature st; 




T(K) 

FIG. 3: (color online)Specific heat vs temperature at constant magnetic fields H. A low temperature A- 
anomaly is evident at i7 > 28 T. The anomaly first moves to higher temperatures with increasing magnetic 
fields, it reaches a maximum at i7 = 36 ± 1 T, and then decreases for 37 < i7 < 45 T. Inset: Specific heat 
vs T at i7 = 37T after subtraction of a small contribution due to the 5^ = triplet level only relevant at 
higher temperatures, and phonons. Also displayed is the result of our Monte Carlo calculation. 

The observation of a phase transition is further confirmed by specific heat measurements as a 
function of temperature for constant magnetic fields shown in Figure 3. This experiment was per- 
formed in a calorimeter made out of plastic materials using both a superconducting magnet 
and the 45 T hybrid magnet at the National High Magnetic Field Laboratory. At zero magnetic 
field the specific heat is featureless. As soon as finite external field along the c-axis is applied 
the specific heat increases due to the reduction in the singlet-triplet energy gap A (77). A small 
anomaly then develops for 77 = 28 T at T ^ 2.7 K that moves to higher temperatures as the 
magnetic field is further increased. This A-shaped anomaly also grows with magnetic field due to 
the increased number of triplons (magnetization). Their number at low temperatures is roughly 
proportional to the magnetic field (see Figure 2(a)). Once the middle of the magnetization ramp is 
reached at 77 = 36 ± 1 T, the transition temperature and the size of the A-anomaly then start to de- 
crease owing to the reduction in the availability of singlet dimers. Moreover, the mirror symmetry 
of the phase diagram around this maximal point (see Fig. 4) offers an experimental confirmation 
of the particle-hole symmetry implicit in 77^^. No phonon contribution has been subtracted from 
the specific heat data. The inset of figure 3 displays the data at 37 T after subtraction of a small 
exponential contribution (activated energy A(77 = 0) = 3.13 meV) from the = Q triplet com- 
ponents that becomes apparent only above 8 K. In addition, a Debye phonon contribution with 
Qd = 350 K was subtracted. The solid line is the result of our quantum Monte Carlo calculation 
for 77^^ after a finite size scaling to the thermodynamic limit. The systems sizes that were used for 
the scaling are: L = 4, 6, 8, 12. 

The phase diagram of figure 4 shows all our experimental data combined, together with some 
temperature traces measured while sweeping the magnetic field at a rate of about 12T/min, i.e. 
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FIG. 4: (color online) Left y-axis: Magnetization normalized to the saturation value (M/Mgat) vs magnetic 
field along the c-axis, measured at 1.5 K (red line). Results for the model discussed in the text (green line). 
Right y-axis: Transition temperature from specific heat vs temperature, magnetocaloric effect (MCE) and 
magnetization vs field data. Black lines are the sample temperature measured while sweeping the magnetic 
field quasi-adiabatically. 

rn 

magneto-caloric effect (MCE) [24]. The left side y-axis is for the magnetization versus H curve. 
From this curve we obtain: Hd = 23.5 T and Hc2 = 49 T. The right side y-axis is for the 
transition temperature versus field plot obtained from the specific curves (see figure 3) and the 
magneto-caloric effect measurements. The anomalies in the MCE, observed upon crossing the 
phase boundaries in the direction of increasing fields, evidence reduced magnetic entropy within 
the critical region consistent with a gas-to-liquid phase transition. Also displayed in figure 4 are 
the transition points calculated with our Monte Carlo algorithm for L = 12. The BEC phase is 
represented by the shaded region whose boundary is intended merely as a guide to the eye. The 
excellent agreement between experiment and theory is a compelling case for BaCuSi206 being a 
true realization of a quasi two dimensional BEC of magnetic degrees of freedom. 

Early Chinese chemists could not have imagined two thousand years ago that BaCuSi206 was 
not only an attractive purple pigment but also a potential solid state device for exploring the quan- 
tum effects of a BEC at liquid ^He temperatures in magnetic fields. Indeed, we show that a minimal 
Hamiltonian for a lattice gas of hard core bosons describes surprisingly well the magnetization and 
specific heat of this magnetic system in external magnetic fields. To derive an effective Hamilto- 
nian we simply neglected the high energy triplet states \(j)l,S^ = 0) and \(j)l,S^ = —1). The 
excellent agreement with the experimental data, and the verification of the particle-hole symme- 
try implicit in Heff validate our approach. However, other physical properties such as the small 
staggered magnetization that appears in the ab plane, when the triplets condense, would require 
the inclusion of virtual processes that involve the above mentioned states. It is also important to 
remark that the phase coherence of the spin BEC requires the conservation of the number of par- 
ticles that is equivalent to the total magnetization. This condition is fulfilled whenever the crystal 
has some degree of rotational symmetry about the direction of the applied field like, for instance. 
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a four-fold symmetry axis. Otherwise, the presence of anisotropy terms will significantly reduce 
the de-coherence time of the condensate. Finally, the weakly coupled Cu-bilayers in the crys- 
tal structure of BaCuSi206 open the door for future studies like the dimensional crossover of a 
bosonic gas. By increasing the distance between adjacent bilayers with chemical substitutions the 
A-shaped second order phase transition, characteristic of a three-dimensional system, is replaced 
by a Kosterlitz-Thouless phase transition. 
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